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ABSTRACT 



We introduce a model with conserved dynamics, where nearest neighbor pairs of 
spins ti (it) can exchange to assume the configuration |t (Ti)? with rate /3 (a), through 
energy decreasing moves only. We report exact solution for the case when one of the 
rates, a or /3, is zero. The irreversibility of such dynamics results in strong dependence 
on the initial conditions. Domain wall arguments suggest that for more general models 
with steady states the dynamical critical exponent for the anisotropic spin exchange is 
different from the isotropic value. 
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1. INTRODUCTION 



Dynamics of one-dimensional systems has received much attention recently. 
One-dimensional models can be solved exactly in some cases, and they provide con- 
venient test cases for theories of reactions, deposition, ordering. Furthermore, in many 
instances, such as, for instance, diffusion-limited reactions, low-dimensional models of- 
fer examples on non-mean-field fiuctuations.^ The scope of models studied recently has 
been extended significantly^ from the "classical" Glauber^ and Kawasaki^ models of the 
Ising-spin flips and exchanges. 

One of the interesting recent developments has been the introduction of anisotropic 
(spatially biased) dynamical moves in irreversible models of reactions. Anisotropic 
particle motion also plays key role in hard-core-particle "asymmetric exclusion" models. 
Steady states of kinetic models with anisotropy have been studied extensively in the field 
of driven diffusive systems^^. Dynamics without detailed balance has also been consid- 
ered recently in the field of neural networks. Our aim in this article has been to initiate 
investigation of the effect of making pair exchanges anisotropic in Kawasaki-type spin 
exchange models of dynamical behavior. 

Most of our results apply in the fully irreversible limit, with only energy-lowering 
moves allowed. This corresponds to the zero-temperature limit and, in the case of 
Kawasaki spin exchanges, yields models of freezing. Studies of irreversible low-temperature 
kinetic Ising models, both spin-flip, Glauber-, and spin-exchange, Kawasaki-type, have 
attracted attention recently.^ The emphasis has been on derivation of exact results for 
these and certain related reaction models. In particular, models with conserved order 
parameter, and related systems, were solved by a method which will be also employed 
here.-*^^"^*^ Various exactly solvable one-dimensional models of freezing have been treated 
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by other approaches. ^-"^ 

Contrary to non-conserved dynamics models of the Glauber type (spin flip), where 
several exact results are available, conserved dynamics, Kawasaki type models have long 
resisted any direct analytical solution. Few exact results have become available recently 
for variants in the zero temperature limit. The ordered domain structure in these 
T = models becomes frozen at large times. 

Irreversible dynamics of the "freezing" type usually retains infinite memory of the 
initial conditions. We will use two different initial conditions, namely, the fully alter- 
nating lattice corresponding to a concentration of ^ of both spin types, and random 
initial conditions where each site is either occupied by a t spin with probability p or by 
a I spin with probability q = 1 — p. 

The outline of the paper is as follows: We start by defining the model in Section 2. 
Then in Section 3 we present the exact solution for the single rate (fully anisotropic) 
case. Finally, in Section 4 we discuss the results and consider generally the anisotropic 
exchange models in the framework of domain-wall arguments which, for the T > 
dynamics, yield information on the dynamical critical exponent. ^""^'^^"^^ 
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2. THE MODEL 



When two spins are exchanged in a ferromagnetic-interaction model, the total en- 
ergy can be increased, unchanged, or decreased. By ferromagnetic interaction we mean 
that the energy of a "bond" connecting parallel spins (tt or J, J,) is lower than the energy 
of a bond connecting antiparallel spins (t| or 

In order to simplify the notation, we will use the language of particle exchanges 
from now on. Thus, the t (i) spins will be replaced by particles A (B). We consider only 
the dynamics in the T — > limit and allow for energy-decreasing particle exchanges. 
Such dynamics leads to frozen states. ■'^'^"■'^^'^^ 

Consider a one-dimensional lattice where each site is occupied by particle A or B. 
The only two possible nearest-neighbor particle exchanges that lower the energy locally 
are 

■■■ABAB--- ^ ■ ■■AABB--- (rate a) (1) 

• • • BAB A ■■■ ^ ■■■ BBAA ■ ■ ■ (rate /3) (2) 

Note that the particle (spin) exchanges are anisotropic. The usual isotropic Kawasaki- 
model spin exchanges correspond to the equal rates a — p. For finite temperatures, 
exchanges other than those shown in (l)-(2) are allowed, and their relative rates depend 
on the temperature.^ In the Ising model notation, exchanges (l)-(2) lower the energy 
by 4J, where — J is the energy of the AA and BB bonds, while J > is the energy of 
the AB and BA bonds. 

The model defined by the moves (l)-(2) is conveniently analyzed^^ in terms of the 
probabilities Pc,n{t) that a randomly chosen interval of n > 1 lattice sites is fully alter- 
nating, ABAB • • • or BAB A • • •, at time t. Here C e {A, B} denotes the starting site 
type in the n-site sequence. For n > 3, the interior of such intervals is "fully reactive," 
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i.e., all internal nearest- neighbor exchanges are allowed (reduce energy). In order to be 
able to claim that the probabilities Pc,nit) are independent of the position along the 
one-dimensional lattice, we must assume that the initial conditions are uniform; specific 
choices will be specified later (see also Section 1). It is important to emphasize that 
the probabilities only refer to the interior of the interval which might or might not be a 
part of a larger alternatively-ordered interval. The probabilities are not conditioned on 
the arrangement of particles outside the selected n-site sequence. 

A closed hierarchy of rate equations can now be written as follows. Consider first 
the case of even n = 2A;, where A; > 1. The rate equation for PA,2k{t) is obtained in the 
form 

-PA,2k = {k- l)aPA,2k + {k- 2)/3PA,2k + f3PA,2k+l + f3PB,2k+l + 2aPA,2k+2 (3) 

Here the dot denotes the time derivative. The first term in (3) represents the rate at 
which the k — 1 sequences of the type (1), internal to the selected interval of 2k consec- 
utive sites, undergo the process given in (1) thereby destroying the fully alternatively- 
ordered state in the selected interval. Similarly, the second term corresponds to the 
k — 2 different four-site subsequences in the selected interval, which undergo process 
(2). 

The third and the fourth terms correspond to, respectively, the rightmost and the 
leftmost AB pairs in the selected interval of 2k sites, exchanging as parts of four-site 
sequences of type (2). In such a process the four-site sequence containing the end-pair 
of the selected interval must in fact be a part of a larger, (2/c-|-l)-site, ordered sequence. 
Thus, Pc,2k+i is used. Finally, the last term corresponds to processes in which the two 
end sites of the selected interval exchange with sites external to it. The configuration 
is of type (1) at both ends, and the alternating order must extend to a larger interval 
of 2/c -|- 2 sites, in order for process (1) to occur. 
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Similarly, we get 



-PB,2k = {k- l)pPB,2k + ik- 2)aPB,2k + aPB,2fc+l + aPA,2fc+l + WPB,2k+2 (4) 

It is also of interest to write down the rate equation for the case n = 2. Indeed, the 
quantities Pc,2 yield densities of AB and BA interfaces in the system. The interfaces 
separate ordered domains of A and B particles (t and [ spins). We have 

-Pa,2 = 2aPA,4 + PPba (5) 

-Pb,2 = WPb,A + OiPA,A (6) 

We next turn to the case of odd n = 2A; + 1, with A; > 1. The rate equations now 

read 

-PA,2k+l = {k- l)aPA,2k+l + {k- l)PPA,2k+l + PPB,2k+2 + OiPA,2k+2 + (« + /3)^'A,2fe+3 

(7) 

-PB,2k+l = (/C - l)/3-PB,2fc+l + {k- l)aPB,2k+l + OlPA,2k+2 + PPB,2k+2 + (« + /3)PB,2k+3 

(8) 

We must complement the above system of rate equations with initial conditions. 
For the initially alternating lattice of A and B particles, the probability of starting at 
either site type is 1/2 regardless of the interval size, i.e., 

<2fc(0) = <2/c(0) = <2fc+l(0) = <2/c+l(0) = I (9) 

For the initially random distribution, placing A particles with probability p and B 
particles with probability q — 1 — p, aX each lattice site, the initial conditions are 

PaUQ) - PBr2kio) - (pq)" (10) 

pr^k+M - p(p^)' (11) 
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The above system of rate equations turns out to be exactly solvable in some cases. 
The case a = (3 (isotropic exchanges) has been considered in the literature. ^^"^^'^"^ We 
could not find a solution for the general case of both rates being non-zero and not 
equal to each other. In the next section we solve the case of a > and /3 = 0. We 
note, however, that the transformation used to eliminate the /c-dependence (see the next 
section) can be generalized to the case of a > 0, /? > 0. Thus, a system of ordinary 
differential equations can be obtained and analyzed numerically or by approximation 
methods. 
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3. THE SINGLE RATE CASE 



In this section we consider a solvable case of the fully anisotropic exchanges. Specif- 
ically, we take /? = 0, a > 0, so that only exchanges of the type (1) are allowed. It 
is furthermore convenient to redefine the time variable (or alternatively, set a = 1) in 
such a way that 

*new = a^old (13) 

This eliminates the ct-dependence from the rate equations. 

More importantly, the /c-dependence can be also eliminated, by a transformation 
which generalizes the one used in the isotropic case.^^~^^'^^ Thus, we try a solution of 
the form 

PAMt) = S{t)^'{t) (14) 
PsMt) = U{t)^\t) (15) 

PA,2k+l{t)^V{t)^''{t) (16) 

PB,2k+i{t)-W{t)^'it) (17) 
Direct substitution in (3)-(4), (7)-(8) yields (with a = 1, /? = 0) 

7 - -7 (18) 

S={1- 2^)S (19) 

if = 2U -{V + W) (20) 

V = (^l-^)V --fS (21) 

I^ = (1 - - -fS (22) 
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We must also check that the proposed forms (14)-(17) are consistent with the initial 
conditions. Denoting the initial values by the subscript 0, we identify 

7o" = 1 (23) 

and 

iT = m (25) 

5™" = U^^"" = 1 (26) 

Vr = P (27) 

WJ^° = q (28) 

These two types of initial conditions were introduced in Section 2. 

For both initial conditions the solution can be now obtained straightforwardly. 
First, we solve (18), 

7(t) = 7oe-* (29) 
This result is then used to solve (19), 

Sit) = So exp [27oe-* - 270 + t] (30) 

The next step is to solve the equations for V and W, (21)-(22), 

Vit) = (Vo - So) exp [7oe-* - 70 + ^] + ^0 exp [2706"* - 270 + t] (31) 

W{t) = (Wo - So) exp [7oe-* - 70 + t] + 5o exp [2706"* - 270 + t] (32) 
Finally, the equation for U is solved, 

U{t)= [Uo e ' H exp 706 - 70 + 2t 

V 70 / 70 ' 
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+ — exp [27oe-* - 270 + 2t] (33) 
70 

These expressions provide the fuU, time-dependent solution of the problem. Specif- 
ically, the probabilities Pc,n{t) are obtained from (30)-(33), with (14)-(17), for all n>3. 
However, for n = 2 special expressions apply, 

PaM^) = ^A,2(0) + 5o7o exp [27oe-* - 270] - SqJo (34) 

PB,2{t) - PbAO) + ^ exp [27oe-* - 270] - ^ (35) 
These results were obtained by integrating (5)-(6), with a = 1, /? = 0. 
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4. DISCUSSION 



As in the isotropic-exchange case, ' the most profound feature of the dy- 
namics is its strong dependence on the initial conditions. In fact, the initial values 
are "remembered" for all times, and enter the t — > oo limiting expressions. The large- 
time limiting expressions can be obtained by expanding in powers of the quantity 7oe~* 
(= 7(t) ) which enters in the exponentials in (31)-(35). In fact, for large times we have 
S,V,W (X e*, U oc e^*, where the proportionality constants depend on the initial values. 
Thus, only the shortest alternatively-ordered intervals survive at large times. Specifi- 
cally, Pa,2(oo), Pb,2(oo), Pa,3{oo), Pb,3(oo), Pb,4(oo) "survive" (remain nonzero) for 
the fully anisotropic dynamics. All other probabilities vanish exponentially. 

The final configuration is therefore incompletely ordered. Some interfaces remain at 
t = oo. These are pairs AB and BA, the probability of which is given by Pa,2 and Pb,2, 
respectively. Some of the interfaces survive as nearby pairs, represented by the nonzero 
values of Pa, 3 and Pb,3- Furthermore, some triplets of interfaces remain unreacted, due 
to /? = 0, corresponding to configuration (2). The probability of the latter is given by 
-Pb,4- For illustration, let us calculate the total density of interfaces per site, I{t), in 
the system, which was also obtained in the isotropic-model studies. We get 

I{t) = PA,2{t) + PbA^) (36) 
/ait(^)= l|l + 3exp[2(e-^-l)]} (37) 
^ran(^) = ^ { 1 + 3 exp [2pq (e-* - l)] } (38) 

We note that both the large-time values themselves, and the coefficients of the expo- 
nential approach to these values, depend on the initial conditions. 



For the initially fully alternating lattice, it is interesting to notice that 

PlMt) = PX%t) = P^S) = Pi^it) = \ exp [2 (e-* - 1)] (39) 

For this particular case, the lattice can be viewed as initially fully covered by AB 
"objects." Each nearest-neighbor pair of such objects can "react away" via the process 
(1). At the end of the process, all the original unreacted AB pairs remain in the centers 
of configurations BAB A. This process is therefore identical to Random Sequential 
Adsorption of dimers on the linear lattice,^® and in fact one can derive the expression for 
PB,'i{t) from the exact expression available for the adsorption of dimers. However, other 
quantities in (39) are less straightforward; we were not able to propose a "dynamical" 
argument for the equalities (39). 

We turn next to domain wall arguments as applied to more general one-dimensional 
models with anisotropic spin exchanges. We will not specify the model but only assume 
that the | spins favor to move to the left, while the J, spins are favorably moved to the 
right. In the particle notation, BA AB exchanges are favored over AB —>■ BA. 

For very low but finite temperatures, the configuration mainly consists of single 
isolated interfaces, separating large ordered domains, at the distance of the correla- 
tion length, ^. If we denote by AE = 2J > the energy cost to create an interface 
by "breaking" a ferromagnetic bond, then for equilibrium states entropy arguments^^ 
suggest that 

C ~ exp [AE/ksT] = exp [2J/A;bT] (40) 

for low temperatures. Of course, with anisotropic spin exchanges, the system may 
not have an equilibrium state but rather reach a steady state^^ at large times. The 
identification (40) then must be checked by exact or approximate methods. 

The usual domain wall argument^^'^^"^^ for isotropic conserved order parameter 
models at low temperatures, can be phrased as follows. The slowest relaxation mode 



in the system will be a locally correlated, difFusional motion of isolated interfaces, on 
the time scale proportional to the inverse of the rate of diffusion and to the separation 
squared, i.e., to However, the rate of diffusion is rather small. An isolated interface 
cannot move on its own by spin exchanges. It must first be turned into a triplet of inter- 
faces by a spin exchange of energy cost 4J. Thus, the diffusion rate will be proportional 
to the inverse Boltzmann factor, exp [— 4J/A;b7']- 

In addition, a pair of interfaces will "unbind" and diffuse away at the rate of order 
1 (because spin exchanges that lead to motion of a pair of interfaces do not cost energy). 
This pair must reach a neighboring isolated interface at a distance of order ^ for the 
process to be complete. If it instead returns to the interface from which it originated, 
then the configuration will be restored to the original one. The probability of not 
returning back is ~ this yields another factor of ^ in the relaxation time.^-*^'^^ 

The resulting expression for the relaxation time r is written as 



where we used (40). This yields the value ^ = 5 for the dynamical critical exponent in 

T OC 

For anisotropic exchange, the unbinding interface pairs will be immediately re- 
turned back for "attracting" interfaces, but they will move away by biased diffusion 
(with asymptotically negligible probability of return) for "repelling" interfaces. (The 
two types of interfaces alternate in the system.) Therefore, the probability factor of not 
returning back, the inverse of 1/^, will not be present. We expect 



'isotropic 



(41) 



Tanisotropic ~ exp [AJ/ksT] 



(42) 



If (40) applies for the specific model's steady state, then we identify 2; = 4, a different 
value as compared to the isotropic case. 



In summary, we reported new exact results for the anisotropic exchange model with 
only the most energy-favorable moves allowed. This restriction allows exact solvability 
in the fully anisotropic case. The resulting dynamics is that of freezing. Domain wall 
arguments suggest that for more realistic, T > models with similar dynamical rules, 
there will be new interesting effects such as the new value of the dynamical critical 
exponent. 
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